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ABSTRACT 



We use a model of an accretion flow coupled with an emergent flare to interpret 
the latest 1.3mm VLBI measurements for Sagittarius A*. The visibility data 
constrained the distances from the flare center to the black hole center as d^w ^ 
20Rg and c^ns ^ 80Rg in the East- West and North-South directions, respectively. 
If interpreted by the hot-spot model, the flare was preferred to pass in front of 
the black hole at a radius much larger than c^ew- If interpreted by the episodic 
jet launched from a nearly edge-on hot accretion flow, the flare was preferred to 
be ejected with 6j > 40° off the black hole rotating axis. This method can be 
generalized to help us understand future sub-millimeter VLBI observations, and 
study the millimeter /sub-millimeter variabilities in the vicinity of the Galactic 
Center supermassive black hole. 



Subject headings: Galaxy: center — submillimeter: general 



1. INTRODUCTION 



Sagittarius A* (Sgr A*) is the be st black hole candidate that can be studied by Very 
Long Baseline Interferometry (VLBI). iDoeleman et al.l (120081 ) reported 1.3mm VLBI detec- 
tions of Sgr A* on two baselines, namely ARO/SMT (Arizona) - CARMA (California) and 
ARO/SMT - JCMT (Hawaii). These data were fitted by a circular Gaussian component 
with a full- width of half maximum (FWHM) of ~ 43/ias, reaching the event-horizon scale 
of the Galactic Center supermassive black hole. Several groups soon interpreted these data 
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Recently, iFish et al.l ( 120111 ) reported new measurements on all the three baselines of the 
ARO/SMT-CARMA-JCMT array over three nights: 2009 April 5-7 (Days 95-97). The data 



on Da y 95 and Day 96 show high consistency with the data in 2007 April by lDoeleman et al. 
( 120081 ) . fitted by a circular Gaussian component with FWHMs of ~ 41/ias and ~ 44/ias 



respectively. But the correlated flux density on all the baselines increased on Day 97. The 
Gaussian model with a similar size but higher flux density is needed to fit these data. 
However, they pointed out that some model parameters must be changed for Day 97 in 
order to fit the data on all the three days by a ring model. This implies that an additional 
variable component might exist together with the accretion flow associated with a black hole 
shadow structure. 

In this Letter, we introduce a method in Section 2 to constrain the position of the 
flaring component, as required to interpret the Day 97 data. In Section 3, we hold general 
discussions on the resolution of the baselines, challenges to the flare models, and suggestions 
for future observations . Summary is pro vided in Section 4. Here we adopt black hole mass 
Mbh = 4.1 X lO^M© (IGhez et al.l boOSi ) and distance of the Galactic Center doc = 8kpc. 
Thus 1 milli-arcsecond (mas) angular size corresponds to a linear size at the Galactic Center 
of 200Rg, where Rg 



GMbh/c^ is the gravitational radius. 
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CONSTRAINING THE POSITION OF FLARING COMPONENT 
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Fig. 1. — Left: Non-scattered image of accretion flow associated with a black hole shadow 
structure (Q Model). Right: Corresponding visibility amplitudes Sq (black lines) and closure 
phases 0cl (plus signs) of baselines I, II and III. Visibility data are shown in squares for Day 
95 and in triangles for Day 96. 



Hereafter, for simplicity, we denote the three baselines from short to long, namely 
ARO/SMT-CARMA, CARMA-JCMT, and ARO/S MT-JCMT by b aselines I, II, and III, 
respectively. We apply all the data on Days 95-97 (IFish et al.l 1201 ll ) and take averages if 
there are more than one data points measured on the same baseline coordinates {u, v) on 
the UV plane. We interpret the measurements on Day 95 and Day 96 as contributions by an 
accretion flow in its quiesc ent state. The mag neto-rotational-instability (MRI) dominated 
relativistic accretion flow ( iHuang et al.l l2009b( ). named as Q Model, with black hole spin 
a = 0.5, inclination angle i = 60°, and position angle = 100°, is adopted to account 
for the emission in the quiescent state. The black hole spin and orientation are typical 
i n our earlier work, which are also in agreement with the es timates made by other groups 
toexter et al.lboiol : Ishcherbakovlboid ISroderick et al.lboill ). 



We simulated the image (shown in the left pa nel of FigJT]) of the accretion flow by ray- 



tracin g with polarized relativistic radiative transfer (iShcherbakov fc Huangll2011t iHuang fc Shcherbakov 

20111 ). We further plot the visibility d ata of this theor etical image in right panel of FigJH 

with inter-stellar scattering considered (IShen et al.ll2005l ). The total flux density Sqo = 2.1Jy. 

The amplitude of visibility, i.e. the correlated flux density Sc shown in black lines flt the 

data on Day 95 with reduced chi-squares Xdof = 1-06 and those on Day 96 with Xdo{ — 1-07, 

here dof = 12 on both days. The closure phases 0cl of baselines I, II, and III averaged over 

10 minutes, as shown in plus signs, is in the range (0°, 25°), which is also consistent with the 
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observations at 0° ± 40° flFish et al.lboilh 
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Fig. 2. — Left Panels: Non-scattered images of accretion flow coupled with a flare with 
different central shifts. Boundaries of preferred positions of flare are shown in dashed white 
lines. Right Panels: Visibility amplitudes and closure phases corresponding to left images. 
Absolutes of 0cl are taken for a concise layout. Visibility data are shown in circles for Day 
97. Cases shown in the top and bottom panels are consistent with the observations, while 
case in the middle panel is not. See text for more details. 



We ascribe the increase of visibility on Day 97 relative to Days 95-96 to the emergence of 
a flaring component. For simplicity, we assume that the brightness distribution of this flare 
can be represented by a circular Gaussian function i3p(Spo, Dp, xp, yp) with four parameters, 
namely Spo (total flux density). Dp (FWHM), xp (shift of flare center) in abscissa, and yp 
(shift in ordinate). The brightness distribution of the flaring state of the accretion flow. Bp, 
is simply assumed to be the superposition of the brightness distribution of the quiescent 
state, Bq, and that of the flare. Bp, i.e. Bp = Bq + Bp. The visibility distribution with 
inter-stellar scattering considered, Vp = {Bp^sca.t}, is also the sum of Vq = {Bq^scat} and 
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Vp = {i3p,scat}, according to the property of linearity of Fourier Transform T. By the 
property of translation of Fourier Transform, we further obtain 

Vf = Vq + Vpo 

Vpo = ^{Sp,scat(Spo,Dp,0,0)}e-2-(-p-^+^p-^), (1) 
where u and v are the baseline coordinates on the UV plane. 

We make a search over wide ranges of the parameters, Spo, Dp, xp, and ?/p, and fit the 
data on Day 97 to calculate the reduced chi-square Xdof with dof = Ai'data — A^param = 14 — 4 = 
10. We choose sets of parameters that yield Xdof = 1 ± 1 by preference. Thus, we obtain 
Spo ~ 0.95±0.55Jy and Dp ~ 6.2±6.2Rg. We estimate the observational angular size of the 
flare as 6'obs = SSl^^mas with the formula = (Dp ■ + 9'^^^^, where ^Rg ~ 0.005mas is 
the angular size of gravitational radius an d ^srat ^ 0.021mas is the extrapolated angular size 



of scattering screen in the EW direction ( IShen et al.l 120051 ). The corresponding brightness 
temperature of the flare is = 1.22 x 10i2(^So/Jy)(z//GHz)"^(^obs/mas)"^ ~ l.Stg;^ x 10^°K. 
This value is consistent with the estimation for a pl asmon adopted to interpret the variability 



of Sgr A* on 3mm wavelength observed by ATCA (ILi et al.l 120091 ). We show the boundaries 
of the preferred positions in dashed white lines in the left panels of Figj2l i.e. all the 
preferred pairs of {xp,yp) located in the region surrounded by the boundaries. This region 
is narrow along the direction of ~ 65° east-of-north, hereafter denoted as EW direction, 
and elongated along the direction of ~ 25° west-of-north, hereafter as NS direction. The 
region has an asymmetrical shape that slightly varies with the parameters of the Q model. 
Generally speaking, the preferred distances from the flare center to the black hole center are 
^Ew ^ 20Rg and c^ns ^ 80Rg, along the EW and NS directions, respectively. 

The data on Day 97 impose a strict constraint on (igw but a relatively relaxed one on 
d^s- This is mainly because the measurements on the tracks of the two long baselines II and 
III are roughly aligned in the EW direction. These baselines are longer than 3GA, giving 
a high angular resolution of 85/ias, or 17Rg, in the EW direction. A flaring component 
emerging at c^ew ^ l^Rg would cause a null/valley point in the visibility profile at a baseline 
length < 3GA, so that the amplitudes of visibility on baseline II might be lower than those 
on baseline III. However, the amplitudes actually have a monotonically decreasing profile 
from baseline II to baseline III, which implies there is no structural variability with scales 
larger than 17Rg in the EW direction. Moreover, the track of baseline I can also detect 
the separation of the two components in the EW direction sensitively, although its precise 
direction is ~ 70° west-of-north, which is ~ 45° off those of baselines II and III. On the 
contrary, d^s cannot be well-constrained since there are no long enough projected baselines 
available in the NS direction. 



In the left panels of Figj2l we show three examples of the flaring state, each consisting 
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of the accretion flow overlaid by a different flare. Pi, shown in the top-left panel, with 
Spio = 0.8Jy, Dpi = 6.6Rg, and (xpi, ypi)/Rg = (—4.08,-2.04), is located within the 
preferred region. The visibility corresponding to Q + Pi is shown in the top-right panel, 
fitting the data on Day 97 by Xdof ~ 1.13. The closure phase (pch is within the range 
(—20°, 20°), indicating high symmetry of the total image. P2, shown in the middle-left 
panel, with Sp2o = O.SJy, Dp2 = 6.6Rg, and (a;p2, |/P2)/Rg = (—12.24, —6.12), is outside the 
preferred region. The visibihty fits the data by Xdo( — As shown in the middle-right 

panel, a flare slightly beyond the boundaries in the EW direction can decrease Sc a lot on 
baseline I and cause a valley structure on baseline II. 0cl is in the range (—40°, 40°), deviating 
greatly from zero, and inferring that the symmetry of the total image in the EW direction 
is somehow broken. P3, shown in the bottom-left panel, with Sp3o = 0.95Jy, Dp3 = 6.2Rg, 
and (iCpa, yp3)/Rg = (—16.34,34.68), is within the preferred region. This is the case with 
minimal Xdof ~ 0-25, i.e. the best-fit. As shown in the bottom-right panel, it even reproduces 
the ~ 10% increase in Sc related to baseline length on baseline I. However, considering the 
uncertainties of the data and the simplicity of the static models adopted here, we think it 
may have over-fitted the data, especially on baseline I. 
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3. DISCUSSION 

We introduce a simple but u seful method to understand the time- variable emission of 



Sgr A* detected by 1.3mm VLBI (IFish et al.ll201ll ). We interpret the data on Days 95-96 by 



an accretion flow in its quiescent state, and the increase of correlated flux density on all the 
three baselines on Day 97 by the emergence of an extra flare, modeled by a circular Gaussian 
spot. The visibility measurements impose a strict constraint on the distance from the flare 
center to the black hole center in the EW direction (65° east-of-north) as c^ew ^ 20Rg. 
They also place another relatively relaxed constraint on the distance in the NS direction 
(25° west-of-north) as c^ns ^ 80Rg. In the rest of this Section, general discussions on various 
aspects will follow. 



3.1. Resolution of Baselines 

The longest baseline length on the tracks of baseline III (ARO/SMT-JCMT) is ~ 
3.55GA, corresponding to the highest resolution of 6'bcam ~ 70/ias. The apparent size of 
the flaring component assumed in this Letter is < 66/x as, slightly super- resolved by the cur- 



rent baseline. Here, we follow the criterion provided in lShen et al.l (119971 ) to judge the degree 
of resolution. 

^LiM, the hmit of source size that can be resolved, is defined as 

^LiM = [6'LiM(Statistical)'^ + ^LiM(Systematic)'^] ^^'^ , 

where 

0.53 

6'lim (Statistical) = --^=== ■ 6'beam 

0.53 

6'lim (Systematic) = —;====-■ 6'beam- (2) 

\/\Fp//\F^\ 



With a mean signal-to-noise ratio (SNR) of 5.5 on the ARO/SMT-JCMT baseline ( iDoeleman et al 



20081 ) and \^Fy/Fy\ being the sum of the uncertainties of the flux density measurements 
~ 30%, this limit is estimated to be 6'lim ~ 20/ias. Thus, we think the flare was resolved in 
an optimistic manner. 
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3.2. Challenges To The Flare Models 

The origin of flares in Sgr A* is still controversial. Various models can be put into two 
general classes, namely the whole change and the transient structure. The whole change 
means enhancement in flux density in the whole emission region, caused by changes in 
physical quantities of the accretion flow, e.g. increase in heating coefficient or creation o f 



power-law electrons (e.g. lYuan. Quataert. fc NarayanI l2003l : iLiu. Petrosian. fc Melial l2004j ). 



We can fit the Day 97 data by the Q Model itself with xLi = 1-07, dof = A^'data = 14, if 
the total flux density is scaled up to 3Jy. This implies those models of whole change are 
plausible for the Day 97 data, i.e. no extra component is required for data interpretation. 
This is consistent with the result of c^ew ^ 20Rg, indicating that the flare happened very 
close to, or even inside the emission region of accretion flow. 

We are also interested in challenging models of the transient structure, which means 
an extra component emerged to contribute to structural variability. Generally speaking, the 
Day 97 data do not preclude any flare model if the transient structure is located within the 
preferred flaring region. 

A popular model used to interpret short-time variability of Sgr A* in millimeter and 



near- i nfrared bands is a ho t-spot orbiting in Keplerian angular velocity (e.g. lBroderick fc Loeb 



20051 : iDovciak et al.l l2008l ) . S uch a hot-spot varying in hourly timescale is predicted to be 



detectable by sub-mm VLBI (IDoeleman et al.ll2009l ). We assume that the flare on Day 97 



was caused by a hot-spot orbiting at ry, the projective radius from the black hole center. 
The normal direction of the orbital plane coincides with the projective rotating axis of the 
black hole. According to earlier work, the position angle of Sgr A* is preferred to be in 
the range (90°, 180°) or (-90°, 0°), rather than in the range (0°,90°) or (-180°, -90°). We 
then obtain ry < 40Rg, as calculated from the constraints on the flaring region shown in the 
above Section. The period of the Keplerian orbit is constrained as TKep(r = ry) < 8hr, which 
is insensitive to black hole spin. The observational duration Tobs was ~ 2hr on Day 97, i.e., 
Tobs > 0.25-TKep(r = ry), implying that an apparent light-curve should be detected with a 
nearly edge-on disk assumed. However, this is inconsistent with the observed sustaining high 
flux density. Therefore, we would either exclude this model, or explain this inconsistency 
by two arguments. One is that the measurements of the total flux density had rather large 
errors. The other is that the radius was shortened by projection effect, i.e. the hot-spot 
might be passing in front of the black hole at a radius r^ > 40Rg with a component aligned 
with the line-of-sight. 

We further generalize this model into a hot-spot moving with sub-Keplerian angular 
velocity and radial velocity comparable to the light speed, the same as the velocity of the 
background transonic accretion flow. In this case, the hot-spot would quickly fall into the 
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event horizo n of the black hole and disappear. We integrate the curve of radial veloc- 
ity shown in iHuang et al.l (l2009b| ) to calculate typical timescale of falling at specific start- 
ing radii ry < 40Rg. The lifetime of the hot-spot is constrained as Tfaii(r = ry) < Ihr, i.e. 
Tobs > 2-Tfaii. Similar to the Keplerian-rotating model, we would either exclude this model 
or interpret it as a special case in which the hot-spot started at a larger radius in front of 
the black hole with r^ > 60Rg, and plunged into the event horizon. 



Furthermore, we consider a different model of episodic jets proposed by lYuan et al. 



(120091 ). in which the flare is contributed by an episodic ejection of plasmoid from a hot 
accretion flow. As they calculated, a plasmoid with initial location of ~ lORg to Sgr A* 
can accelerate from rest to ~ 0.8c in 35 min. This predicts that during the 2h observation, 
the plasmoid can travel ~ 250Rg from the black hole center. To take into account both the 
constraint on c/ns and the projection effect by a disk with i > 60°, the plasmoid must be 
ejected in a direction with 6^ > 40° off the rotating axis of the black hole. 

We wish to make a comment on a model of tidal disruption of asteroid explored by 
Zubovas. Nayakshin. fc MarkofB ( 1201 ll ). in which the flare is related to asteroids rather than 
properties of hot accretion flow. They predicted a small size for a flaring region < lAU 
or < 25Rg, which is roughly included in our constraints with boundaries of d^w and c/ns- 
However, the timescale of an asteroid in a parabolic orbit around Sgr A* is estimated to be 
< 1.5hr, which is unfavored by the Day 97 flare. 
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3.3. Visibility Prediction For An Array of Five Stations 
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Fig. 3. — Left Panels: Non-scattered images of accretion flow coupled with a flare with 
different central shifts. Boundaries of preferred positions of the flare are shown in dashed 
white lines. Right Panels: Visibility amplitudes and closure phases (APEX-ARO/SMT- 
LMT, in plus signs) corresponding to the left images. See text for more details. 

The flaring region of Sgr A* can be constrained better if suitable baselines in the NS 
direction are available. Theoretically, the separation of AS* ~ c/ns ^ 80Rg between two 
components, the quiescent accretion flow and the flare, will cause first null/valley point in 
the visibility profile at Lq > 0.25GA in the NS direction. We include two additional stations, 
APEX in Chile and LMT in Mexico, which are promising for contributing to sub-mm VLBI 
observations in the near future. The tracks of seven additional baselines, from short to long, 
ARO/SMT-LMT, CARMA-LMT, APEX-LMT, JCMT-LMT, APEX-ARO/SMT, APEX- 
CARMA, and APEX-JCMT, cover from 0.5GA to 7GA in length close to the NS direction. 
In an optimistic view, the separation of AS ~ c^ns ^ 40Rg could be detected if all the 
baselines work. In the left panels of FigJHl we choose two examples of flares, P4 and P5, both 
with Sp = 0.95 Jy and Dp = 6.2Rg, but with (xp4, yp4) /Rg = (—2.04, 4.9) and (a;p5, j/ps) /Rg = 
(—16.32, 40.8), respectively. The corresponding visibilities are shown in the right panels. The 
image of the Q Model and P4, with AS ~ 5Rg, causes the first valley point at ~ 4GA on 
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the baseline APEX-LMT. The image of the Q Model and P5, with AS* ^ 40Rg, causes the 
first valley point at ~ 0.5GA on the baseline ARO/SMT-LMT, which appears as a drop in 
amplitude between the baselines ARO/SMT-CARMA and ARO/SMT-LMT. At the same 
time, closure phases by a station group including APEX and LMT, also imply the change 
to symmetry caused by the flare. For example, |0cl| by the group of APEX - ARO/SMT - 
LMT are shown in blue plus signs in the right panels of FigjSl 



4. SUMMARY 

The detection of time- variable emission of Sagittarius A* by 1.3mm VLBI indicates 
instability on event-horizon scales. Modeling the flare by a circular Gaussian spot, the data 
constrained the spot size as Dp ~ 6.2 ± 6.2Rg and deviation to the black hole center in 
the East- West direction (65° east-of-north) as dEw ^ 20Rg, both being comparable to the 
size of the black hole shadow (~ lORg). Alternatively, we can interpret this flaring activity 
as enhancement in whole emission region, and we would prefer such an interpretation to 
models of transient structures that require an extra component. This is mainly due to the 
sustaining high flux density in long observational duration, Tobs ~ 2hr. The duration is 
more than one-fourth of the orbiting period if modeled by a hot-spot in Keplerian rotation 
at a preferable orbiting radius < 40Rg. It is also more than twice the falling timescale if 
modeled by a hot-spot with sub-Keplerian angular velocity and high radial velocity, and 
with a preferable initial radius < 40Rg. A hot-spot located at d-Ew cannot maintain high 
flux density during the observational duration, unless it passed in front of the black hole 
at a radius much larger than c^ew so that the observed distance is shortened by projection 
effect. If interpreted by an episodic jet, the ejecting plasmoid cannot be confined within the 
preferred fiaring region with an assumed edge-on accretion fiow, unless the angle between 
the ejecting direction and the black hole rotating axis is greater than 40°. This method 
of visibility analysis can be generalized for future sub-millimeter VLBI measurements. We 
would have a better understanding of the nature of variability in the black hole vicinity 
of the Galactic Center if c^ns; deviation of the fiare to the black hole center in the North- 
South direction (25° west-of- north) , could be constrained more precisely with new stations 
included. 
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